Anisotropic dielectric function in polar nano-regions of relaxor ferroelectrics 
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The paper suggests to treat the infrared reflectivity spectra of single crystal perovskite relaxors as 
fine-grained ferroelectric ceramics: locally frozen polarization makes the dielectric function strongly 
anisotropic in the phonon frequency range and the random orientation of the polarization at nano- 
scopic scale requires to take into account the inhomogeneous depolarization field. Employing a simple 
effective medium approximation (Bruggeman symmetrical formula) to dielectric function describing 
the polar optic modes as damped harmonic oscillators turns out to be sufficient for reproducing 
all principal features of room temperature reflectivity of PMN. One of the reflectivity bands is 
identified as a geometrical resonance entirely related to the nanoscale polarization inhomogeneity. 
The approach provides a general guide for systematic determination of the polar mode frequencies 
split by the inhomogeneous polarization at nanometer scale. 

PACS numbers: 78.30.-j, 77.80.-e, 64.70.Kb, 77.80.Bh, 



In recent years, there has been enormous effort in 
studying single crystals with intrinsic nanoscopic inho- 
mogeneity, since such materials often show a very in- 
teresting properties. It was even proposed that the 
clustered, inhomogeneous states encountered for exam- 
ple in high-Tc cuprates, CMR manganites, nickelates, 
cobaltites, diluted magnetic semiconductors or ferrolec- 
tric relaxors, should be considered as a new paradigm in 
condensed matter physics X In case of relaxors, the pecu- 
liar dielectric properties of relaxor materials were associ- 
ated with the presence of small polar clusters - so called 
polar nano-regions (PNR's)- already in the pioneer work 
of Burns and Dacol 2 . However, because of their small size 
and random nature, we still lack a clear understanding of 
their size distribution, thickness and roughness of their 
boundaries, their connectivity, shape anisotropy, distri- 
bution of the associated dipolar moments, their fractal 
self-similarity, their dynamics, and so on. PNR's are 
often represented as small islands submerged in a non- 
polar matrix, possibly appearing and disappearing again 
in time. On the other hand, the recent piezoelectric scan- 
ning microscopy investigations^ of the surface of PbTi0 3 - 
doped relaxors rather invoke a picture of a fine, hier- 
archical and essentially static "nanodomain" structure. 
It strongly suggests that the common perovskite relax- 
ors are actually quite densely filled by quasi-static polar 
nano-regions, and that the former picture with prevailing 
non-polar matrix can perhaps be appropriate only at high 
temperatures around the so called Burns temperature^ 

Throughout this paper we will have in mind com- 
mon perovskite relaxors like Pb(Mg 1 / 3 Nb 2 /3)03 (PMN), 
Pb(Zn 1/3 Nb 2/3 )0 3 (PZN), Pb(Sc 1/2 Ta 1/2 )0 3 (PST), 
(Pb,La)(Zr,Ti)0 3 (PLZT) and similar systems. Vari- 
ous experiments show that the dipolar moments in these 
relaxors are caused mainly by ionic off-center displace- 
ments. It is difficult to get reliable information about 
the directional distribution of these displacements, but 
the amplitude of the relevant ion displacements (eg. Pb 
in PMN) is quite well defined^*& and it is of the same 
order as in normal ferroelectrics. We will assume that 



these local displacements are more or less parallel within 
each PNR (one would be in trouble how to define PNR 
if it were not the case) and that the PNR's are at the 
time scale of our interest essentially staticji (i.e. the ions 
vibrate around their displaced but fixed positions, except 
perhaps for those at PNR boundaries) . Under such con- 
ditions, the homogeneous frozen polarization (dipole mo- 
ment density) Pp can be well defined within each PNR, as 
well as the locally homogenous dielectric function e(w), 
describing the contribution of polar vibrations inside a 
given PNR. 

It is obvious that the cubic (and in harmonic approxi- 
mation isotropic) environment of ions in perovskite struc- 
ture is broken by their off-center displacements. Since 
the displacements are aligned within a given PNR, par- 
allel (|| Pp) and perpendicular (_L Pp) ionic fluctuations 
are strongly inequivalcnt. In fact, it is quite likely that 
the parallel fluctuations feel a more stiff potential, as in 
the case of usual ferroelectrics. Within a given PNR, 
one may thus expect that the ionic contributions make 
the e(uj) tensor strongly anisotropic. For simplicity, we 
will assume that the anisotropy is uniaxial, so that the 
e(u>) tensor has only two principal components, eii and e± 
(parallel and perpendicular to Pp). The aim of this paper 
is to demonstrate that infrared (IR) reflectivity spectra 
of common perovskite relaxors can be rather well under- 
stood by taking into account this anisotropy. 

In the view of the above considerations, it seems rea- 
sonable to analyze the influence of PNR on the polar 
phonon modes under the following simplifications: 

• PNR's are frozen at phonon frequencies 

• PNR size is much smaller than the IR wavelength 

• Volume of the sample is fully covered by PNR's 

• Shapes of PNR's are roughly spherical 

• Pp and e(u>) within a given PNR are homogeneous 

• e(w) has uniaxial anisotropy (principal axis || Pp) 
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• Orientations of principal axe are random 

• All PNR's have the same e(w) 

It is known that the reflectivity calculations for the 
samples with dielectric inhomogeneities (spatial fluctu- 
ations of dielectric function) at scales smaller than the 
wavelength of the probing radiation can be performed 
within the effective medium approach^ In this approach, 
the medium is fully characterized by a homogeneous (av- 
eraged) effective dielectric function, e c s, which could be 
in principle evaluated directly from the original heteroge- 
nous dielectric function. In view of the simplifications 
already made, it seems sufficient to use here a rough ap- 
proximation, known as the effective medium approxima- 
tion (Bruggeman symmetrical formula for binary com- 
posites), where e c g is given by the implicit equation 9 



principal bands (below 100, near 220 and 550 cm -1 ) cor- 
respond to the three TO polar modes of the average cubic 
structure. Each of these bands has a tail or a bump on 
the high frequency side. It is natural to assign the prin- 
cipal peaks to the doubly degenerate Ei-components and 
the high frequency wings to the stiffened Ai-components 
of the three polar modes. Similarly, the positions of the 
corresponding three LO bands can be roughly read out 
from the plot of the imaginary part of the inverse per- 
mittivity. 
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The IR reflectivity of a thick (opaque) sample is thcn^ 
evaluated from the standard formula 
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As a matter of fact, we exploit here the equivalence^ of 
our model with that of a dense random binary composite 
of spherical particles with isotropic dielectric functions 
equal to ey and ej_, and relative appearance of 1:2 volume 
ration 

It can be expected that the principal polar modes of 
PNR's in perovskite relaxors are those originating from 
the 3 Fi u polar modes of the parent cubic structure. The 
uniaxial anisotropy induced by the frozen polarization 
split each of these triply degenerate modes in an Ai+Ei 
pair (we assume a strong anisotropy limit leading to Coon 
symmetry for any Pp direction). The principal compo- 
nents of the dielectric tensor can be thus conveniently 
parameterized assuming the factorized formii for gener- 
alized damped harmonic oscillator model: 
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where uj a x , to j and ^Ai.LOj is respectively the transverse 
and the longitudinal frequency of the j-th mode polar- 
ized along the local polarization direction, 7Ai,toj an< i 
7Ai,loj are the corresponding damping parameters and 
e ||,oo is th- e corresponding component of the high fre- 
quency permittivity tensor. The e± tensor component 
determined by Ei modes is defined analogously. 

As an example, let us now apply the model (eqs. (1)- 
(3)) to our room temperature reflectivity data of PMN, 11 
reproduced in Fig. 1. The real and imaginary part of the 
resulting effective permittivity is shown in Fig. 2 and it 
is obviously very similar to that obtained from the stan- 
dard multi-oscillator fit technique^ Clearly, the three 
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FIG. 1: (Color online) IR reflectivity spectra of the PMN 
single crystal. Circles stand for the room temperature data 
of Ref. |Tl] (by mistake denoted as 20K there), full line stands 
for the fit with the model (l)-(3), parameters are given in 
Table I. 

A complete set of mode parameters, as obtained by 
adjusting the model (eqs. (l)-(3)) to the measured re- 
flectivity spectrum, is given in Table 1. The high fre- 
quency dielectric tensor was assumed isotropic, since even 
in PbTiOs its two principal components differ by less 
than 1 %. 12 It was set to eii )0O = ej_ i0 o = 5.75, de- 
duced from the available (average) refractive index data 
of PMN>i^ Furthermore, the visual agreement was much 
improved by adding a weak mode near 335 cm -1 (in- 
troduced as identical in both en and e± ). The agree- 
ment is fairly good and the number of the parameters 
is adequate, except for the frequency of the overdamped 
Ei(TOl) soft mode, which obviously cannot be reliably 
determined from data measured above 20 cm -1 . The 
value selected here makes the total phonon contribution 
to the effective static dielectric permittivity of the order 
of 1000, which is in qualitative agreement with our pre- 
viously published analysis iLl^ilk (There is an additional 
strong relaxation in the GHz frequency region^ which is 
most likely due to the fluctuations of PNR boundaries, 
but it is not in the focus of this paper.) All the param- 
eters appear quite realistic - the mode frequencies are 
reasonably close to that of PbTiOs, given for compari- 
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son in Table 1, and the trends are in agreement with the 
smaller frozen polarization^ in PMN with respect to the 
spontaneous polarization of PbTi03. 
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TABLE I: Intrinsic frequencies and damping parameters (o;,7)~ 
of polar modes of PMN obtained from the fit of the reflectiv- 
ity spectrum shown in Fig. 1 with the model defined by eqs.-r 
(l)-(3) (e||,oo = e±,oo = 5.75), in comparison with the corre- 
sponding mode frequencies of room temperature PbTiU3. 

We do not know for sure the origin of the weak fea- 
ture near 335 cm -1 . In principle, the local homoge- 
neous polarization can induce IR activity also for the 
Ei mode derived from the remaining F2 U (silent) optic 
mode. However, the IR dielectric strength of this mode in 
usual ferroelectric perovskites is known to be extremely 
weak. 12 Another possibility is that this spectral feature is 
related to the random occupancy on (ABO3) perovskite 
B-sites. As a matter of fact, there is also a weak IR ac- 
tive mode expected around this frequencjiii due to the 
B-site short range ordering of NaCl-type, which activates 
the R 25 Brillouin zone corner mode of the parent cubic 
structure (anti-parallel vibrations of incquivalent ions at 
neighboring B-sites). The strength of this mode is deter- 
mined by a compromise between its relatively strong in- 
trinsic strength inside of the ideally ordered regions (Born 
charges at inequivalent B-sites are quite different) and 
relatively small degree of this B-site order. This mode 
was in fact quite clearly seen in IR spectra of compounds 
with a larger degree of B-site order, like PST and similar 
systems with 1:1 B-site stoichiometryjiLiSi so that it is 
perhaps the most likely interpretation. 

The imaginary part of effective permittivity (Fig. 2 
(b)) seems to indicate another weak but clear mode with 
TO frequency near 400 cm' 1 . The origin of this mode 
was a puzzle since realistic calculations^ show that no 
TO mode is expected around this frequency. Interest- 
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FIG. 2: (Color online) Spectra of real (a) and imaginary (b) 
part of the permittivity and of imaginary part of the inverse 
permittivity (c) calculated from the model defined by eqs. 
(l)-(3) with the adjusted parameters given in Table I. Full 
line corresponds to the effective macroscopic permittivity e e g, 
dashed and dotted corresponds to the e± and ey components 
of the local permittivity, respectively. 



ingly, the effective medium model proposed here repro- 
duces this band without assuming any intrinsic TO fre- 
quency near 400 cm . Such modes, so called geometri- 
cal resonances, are actually interfacial modes intimately 
related to the heterogeneity of the medium, and they 
are known for example from IR spectra of ceramics of 
anisotropic materials^ This purely geometric resonance 
related to the double reflectivity minimum near 400 cm -1 
appears to be pronounced due to a relatively large split- 
ting of the L02 mode in comparison with its damping, 
and it is a generic feature present in many other relaxor 
perovkites. 

A systematic analyses of IR reflectivity data of other 
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relaxor perovskites as well as of temperature dependences 
within this model are now in progress. In PMN, the 
frequency of the Ai(TOl) mode increases on tempera- 
ture lowering and it is obvious that this mode coincides 
with that which was denoted soft mode in the previous 
low temperature IR 1114 - 15 and neutron studies w& On the 
contrary, the Ei-component of TOl mode seems to re- 
main rather soft in all perovskite relaxors. Actually, it 
would be extremely helpful to develop an appropriate av- 
eraging scheme also for the polar modes in Raman scat- 
tering and inelastic neutron scattering spectra, in order 
to form a solid basis for comparison. For example, it 
seems to us that Raman bands denoted^ as A-F could 
be related to the polar modes Ai(TOl), Ai(L03), LOl, 
T02 + silent, T03 and L02 modes, respectively. At the 
present stage however, a direct quantitative comparison 
is difficult and could be even misleading. 

The present model can be modified in many ways. For 
example, one may consider an elongation of PNR along 
the direction of the frozen polarization or to consider 
the contribution of the residual non-polar matrix, both 
possible using the analogical effective medium schemes 
developed for ceramicsSS. More interestingly, one can 
hope to incorporate a more realistic information about 
the geometry of PNR structure within a more sophis- 
ticated effective medium approach. At the level of the 
crude approximation (eq. 1), it actually does not mat- 
ter whether the directions of the frozen polarization are 
completely random or whether they tend to be oriented 
preferentially along a family of high symmetry directions 
(most likely, 111 or 100). None of these modifications 
is expected to have a drastic influence on the intrinsic 



phonon frequencies as determined here, but it is possible 
that it may explain the spectra assuming smaller damp- 
ing parameters. 

In conclusion, IR spectroscopy brings evidences of 
splitting of polar optic modes in perovskite relaxors. We 
deduce that at phonon frequencies and at the length-scale 
of PNR's, the dielectric function is strongly anisotropic. 
It is demonstrated on the PMN crystal case that a sim- 
ple effective medium model can reproduce the average 
permittivity and reflectivity spectra fairly well. The ef- 
fective medium approach also explains the pseudo- TO 
mode near 400 cm^ 1 as a consequence of nanoscopic het- 
erogeneity. It is believed that the model can be used for 
systematic estimating the frequencies of the three prin- 
cipal polar modes in perovskite relaxors and their Ai-Ei 
splitting due to the underlying polarization, provided the 
reflectivity data are known in a sufficiently broad interval 
(20-1000 cm" 1 ). 
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